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Magnetic field variations are detected before rupture in the form of 'spikes' of alternating sign. 
The distinction of these 'spikes' from random noise is of major practical importance, since it is easier 
to conduct magnetic field measurements than electric field ones. Applying detrended fluctuation 
analysis (DFA), these 'spikes' look to be random at short time-lags. On the other hand, long range 
correlations prevail at time-lags larger than the average time interval between consecutive 'spikes' 
with a scaling exponent a around 0.9. In addition, DFA is applied to recent preseismic electric field 
variations of long duration (several hours to a couple of days) and reveals a scale invariant feature 
with an exponent a m 1 over all scales available (around five orders of magnitude). 
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Many physical and biological complex systems 
exhibit scale- invariant features characterized by 
long-range power-law correlations, which are of- 
ten difficult to quantify due to the presence of er- 
ratic fluctuations, heterogeneity and nonstation- 
arity embedded in the emitted signals. Here, we 
focus on different types of nonstationarities such 
as random spikes and pseudo-sinusoidal trends, 
that may affect the long-range correlation prop- 
erties of signals that precede rupture. Since these 
nonstationarities may either be epiphenomena of 
external conditions or may arise from the intrinsic 
dynamics of the system, it is crucial to distinguish 
their origin. This is attempted in the present pa- 
per for both the magnetic and the electric field 
variations that appear before rupture by employ- 
ing the detrended fluctuation analysis (DFA) as 
a scaling method to quantify long-range temporal 
correlations. In particular, for the magnetic field 
variations which have usually the form of 'spikes' 
of alternating sign, we find that at short time 
scales they look to be random (thus may then 
be confused with random noise), but at larger 
scales long-range correlations prevail. As for the 
electric field variations, that are usually superim- 
posed on a pseudo-sinusoidal background, upon 
using the longest time series available to date (i.e., 
with duration up to a couple of days), a scale- 
invariant feature over five orders of magnitude 
emerges with an exponent close to unity. 



I. INTRODUCTION 

The detrended fluctuation analysis^— (DFA) is a 
novel method that has been developed to address the 
problem of accurately quantifying long range corre- 
lations in non-stationary fluctuating signals. It has 
been already applied to a multitude of cases including 



diac dynamics^—, meteorolog y 21 ' 22 , climate tempera- 
ture fluctuations 23 . Traditional methods such as power 
spectrum and autocorrelation analysis 2 ^ are not suitable 
for non-stationary signals^. 

DFA is, in short, a modified root-mean-square (rms) 
analysis of a random walk and consists of the following 
steps: Starting with a signal u(i), where i = 1, 2, . . . , N, 
and N is the length of the signal, the first step is to 
integrate u(i) and obtain 



V(i) = 

where u stands for the mean 
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We then divide the profile y(i) into boxes of equal length 
n. In each box, we fit y(i) using a polynomial function 
y n (i) which represents the local trend in that box. (If 
a different order I of polynomial fit is used, we have a 
different order DFA-/, for example DFA-1 if 1=1, DFA- 
2 if 1=2, etc.) Next, the profile y(i) is detrended by 
substracting the local trend y n (i) in each box of length 



Y n (i) = y(i) - y n {i)- 



(3) 



Finally, the rms fluctuation for the integrated and de- 
trended signal is calculated 



F(n) 
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The behavior of F(n) over a broad number of scales is 
obtained by repeating the aforementioned calculation of 
F(n) for varied box length n. For scale invariant signals, 
we find: 
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where a is the scaling exponent. If a = 0.5, the signal is 
uncorrelated (white noise), while if a > 0.5 the signal is 
correlated. 

By employing the DFA method it was foun d 25 ' 26 that 
long-range correlations exist in the original time-series 
of the so called seismic electric signals (SES) activities. 
These are low frequency (< 1Hz) electric signals preced- 
ing earthquakes 2 ^— , the generation of which could be 
understood in the following context. A change of pres- 
sure affects the thermodynamic parameters for the for- 
mation, migration or activation, in general, of defects in 
solids 3 —. In an ionic solid doped with aliovalent impuri- 
ties a number of extrinsic defects is produced ) 34 ' 35 due to 
charge compensation, a portion of which is attracted by 
nearby aliovalent impurities, thus forming electric dipoles 
that can change their orientation in space through a de- 
fect motio n 36 ' 37 . Hence, pressure variations may affect 
the thermodynamic parameters of this motion, thus re- 
sulting in a decrease or increase^ 8 , of the relaxation time 
of these dipoles, i.e., their (re)orientation is taking place 
faster or slower when an external electric field is applied. 
When the pressure, or the stress in general, reaches a 
critical valued a cooperative orientation of these electric 
dipoles occurs, which results in the emission of a transient 
electric signal. This may happen in the focal region of 
a (future) earthquake since it is generally accepted that 
the stress gradually changes there before rupture. 

It has been show n 40 ' 41 that SES activities are better 
distinguished from electric signals emitted from man- 
made sources, if DFA is applied to a signal after it 
has been analyzed in a newly introduced time domain, 
termed natural time %. In a time series comprising N 
events, the natural time Xk = k/N serves as an index^ 5 - 
for the occurrence of the fc-th event. The evolution of 
the pair (xk,Qk) is studie d 25 ' 31 ' 40 " — , where Qk denotes 
a quantity proportional to the energy released in the fc-th 
event. For dichotomous signals, which is frequently the 
case of SES activities, the quantity Qk stands for the du- 
ration of the fc-th pulse. The normalized power spectrum 
n(u>) = |<&(w)| 2 was introduced, where 
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<&(w) = Pk exp 
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and pk = Qk/ Yln=i w = 2tt0; 4> stands for the natu- 
ral frequency. In natural time analysis, the properties of 
n(w) or Tl(4>) are studied for natural frequencies cj) less 
than 0.5, since in this range of <f), H(u>) or H((j>) reduces to 
a characteristic function for the probability distribution 
Pk in the context of probability theory. When the system 
enters the critical stage, the following relation hold s 25 ' 31 : 
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For uj -> 0, Eq.© leads to^ 5 . 

n(w) « 1 - 0.07a; 2 
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which reflects^ that the variance of x is given by 
Ki = (X 2 ) - {x? = 0.07, 



where (/(x)> = J2k=iPkf(Xk)- The entropy S in the 
natural time-domain is defined as 4 ^ 

S= (xln%) - (x) Mx)> 

which depends on the sequential order of events^ 2 -. It 
exhibits^ concavity, positivity, Lesch o 45 ' 46 stability, and 
for SES activities (critical dynamics) its value is smaller— 
than the value S u (= ln2/2-l/4 0.0966) of a "uniform" 



(u) distribution (as defined in Refs4°-~— , e.g. when all 
Pk are equal or Qk are positive independent and identi- 
cally distributed random variables of finite variance. In 
this case, K\ and S are designated n u (= 1/12) and S u , 
respectively). Thus, S < S u . The same holds for the 
value of the entropy obtaine d 43 ' 44 upon considering the 
time reversal T, i.e., Tpk = PN-k+i, which is labelled by 
S-. In summary, the SES activities, in contrast to the 
signals produced by man-made electrical sources, when 
analyzed in natural time exhibit infinitely ranged tempo- 
ral correlation o 40 ' 41 and obey the conditions^ 4 -: 



and 



ki = 0.07 



< S u . 
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For major earthquakes, i.e., with magnitude Mw6.5 or 
larger, SES activities are accompanied 4 ^ by detectable 
variations of the magnetic field B 4 ^. These variations, 
which are usually measured by coil magnetometers (see 
below), have the form of 'spikes' of alternating sign. It is 
therefore of interest to investigate whether these 'spikes' 
exhibit long range temporal correlations. This investiga- 
tion, which is of major importance since only magnetic 
field data are usually available in most countrie s 28 ' 49 ' 50 
(since it is easier to conduct magnetic field measurements 
than electric field ones) , is made here in Section II. 

In the up to date applications of DFA, long-range cor- 
relations have been revealed in SES activities of duration 
up to a few hour s 25 ' 26 ' 40 ' 41 . During the last few years, 
however, experimental results related to some SES ac- 
tivities of appreciably longer duration, i.e., from several 
hours to a couple of days, have been collected. These 
data now enable the investigation of scaling in a wider 
range of scales than hitherto known. This provides an 
additional scope of the present paper and is carried out 
in Section III. A discussion of the results concerning the 
magnetic and electric data follows in Section IV. Finally, 
our conclusions are presented in Section V. 



II. MAGNETIC FIELD VARIATIONS 
PRECEDING RUPTURE 

The measurements have been carried out by three 
DANSK coil magnetometers (DMM) oriented along the 
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three axes: EW, NS and vertical. Details on the calibra- 
tion of these magnetometers can be found in the Supple- 
mentary Information of Ref.—. In particular, this cali- 
bration showed that for periods larger than around half a 
second, the magnetometers measure 5 - 1 - the time derivative 
dB/dt of the magnetic field and their output is 'neutral- 
ized' approximately 200ms after the 'arrival' of a Heavi- 
side unit step magnetic variation. It alternatively means 
that a signal recorded by these magnetometers should 
correspond to a magnetic variation that has 'arrived' at 
the sensor less than 200ms before the recording. The 
data were collected by a Campbell 2 IX datalogger with 
sampling frequency f exp — lsample/s. 

Figure QJa) provides an example of simultaneous 
recordings on April 18, 1995 at a station located close 
to Ioannina (IOA) city in northwestern Greece. Vari- 
ations of both the electric (E) and magnetic field are 
shown. They were followed by a magnitude Mw6.6 earth- 
quake (according to the Centroid Moment Tensor solu- 
tions reported by the United States Geological Survey) 
with an epicenter at 40.2°N21.7°E that occurred almost 
three weeks later, i.e., on May 13, 1995. The record- 
ings of the two horizontal magnetometers oriented along 
the EW- and NS-directions labelled Bew an d Bns, re- 
spectively, are shown in the lower two channels. They 
consist of a series of 'spikes' of alternating sign as it be- 
comes evident in a lOmin excerpt of these recordings de- 
picted in FigJTJb). These 'spikes' are superimposed on 
a background which exhibits almost pseudo-sinusoidal 
variations of duration much larger than Is that are in- 
duced by frequent small variations of the Earth's mag- 
netic field termed magnetoteiluric variations (MT). In 
addition, the horizontal ^—variations were monitored 
at the same station by measuring the variation AV of 
the potential difference between (pairs of) electrodes - 
measuring dipoles- that are grounded at depths of « 2m. 
Several such dipoles were deployed along the EW- and 
NS-directions with lengths (L) a few to several tens of 
meters (short-dipoles) or a couple of km (long-dipoles) 
(Thus the electric field is given by E = AV/L and is 
usually measured in mV/km). For example, the following 
measuring dipoles are shown in the upper three channels 
of Figflja): Two short electric dipoles at site 'c' of IOA 
station (see the Supplementary Information of Ref4£ as 
well as Reb£2 where the selection of site 'c' has been 
discussed) of length 50m (labelled E c -W c and N c -S c , re- 
spectively) while the uppermost channel corresponds to 
a long dipole (labelled L' s -I) with length k 5km at an 
almost NS-direction. As it becomes obvious in Fig[TJb), 
the E variations consist of a series of almost rectangular 
pulses (cf. the initiation and cessation of each rectangu- 
lar pulse correspond to two 'spikes' with opposite sign in 
the B recordings). 

We now apply DFA to the original time-series of the 
magnetic field variations and focus our attention on the 
Bew component where the intensity of 'spikes' is higher. 
Dividing the time series of length N into N/n non- 
overlapping fragments, each of n observations, and deter- 



mining the local trend of the sub-series we find the cor- 
responding log F(n) versus log(n) plot where n — f exp At 
shown in Fig[2J If we fit the data with two straight lines 
(which are also depicted in Figj5J), the corresponding scal- 
ing exponents are a — 0.52 ±0.04 and a = 0.89 ±0.03 for 
the short- and long- time lags (i.e., smaller than around 
12s and larger than « 12s), respectively. The crossover 
occurs at a value (At ~ 12s) which is roughly equal to 
the average duration (T) « 11.01 ± 0.03s of each electric 
pulse, corresponding to the interval between two consecu- 
tive alternating 'spikes'. Thus, FigJ^] shows that at time- 
lags At larger than (T) long-range power law correlations 
prevail (since a > 0.5), while at shorter time-lags the a 
value is very close to that (a = 0.5) of an uncorrelated 
signal (white noise). 

The above findings are strikingly reminiscent of the 
case of signals with superposed random spikes studied 
by Chen et alii. They reported that for a correlated 
signal with spikes, they found a crossover from uncorre- 
lated behavior at small scales to correlated behavior at 
large scales with an exponent close to the exponent of 
the original stationary signal. They also investigated the 
characteristics of the scaling of the spikes only and found 
that the scaling behavior of the signal with random spikes 
is a superposition of the scaling of the signal and the 
scaling of the spikes. The case studied by Chen et al^, 
however, is different from the present case, because the 
'spikes' studied here correspond to the pre-seismic mag- 
netic field variations and hence are not random (cf. recall 
that when applying DFA to the 'durations' of the electric 
field rectangular pulses shown in Figfljb), we found 4 ^ an 
exponent around unity). 



III. DFA OF SES ACTIVITIES OF LONG 
DURATION 

In Fig[3]the following four long duration SES activities 
are depicted all of which have been recorded with sam- 
pling frequency f exp = lsample/s at a station close to 
Pirgos (PIR) city located in western Greece (In this sta- 
tion only electric field variations are continuously moni- 
tored with a multitude of measuring dipoles the deploy- 
ment of which has been described in detail in the Sup- 
plementary Information of Ref. 53 ). First, the SES ac- 
tivity on September 17, 2005 with duration of several 
hours that preceded the Mw6.7 earthquake with an epi- 
center at 36.3°N23.3°E on January 8, 2006. Second, 
the SES activity that lasted from January 21 until Jan- 
uary 26, 2008 and preceded the Mw6.9 earthquake at 
36.5°N21.8°E on February 14, 2008. Third, the SES 
activity during the period from February 29 to March 
2, 2008 that was followed 3 -^ by a Mw6.4 earthquake at 
38.0°N21.5°E on June 8, 2008 (SES activity information 
is issued only when the expected magnitude is around 6 
units or large r 44 ' 53 ). Finally, FigJ^d) depicts the most 
recent SES activity of duration of several hours that was 
recorded on December 12, 2008. The latter SES activity 
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was followed by a Mw5.4 earthquake at 37.1°N20.8°E 
on February 16, 2009 (according to Athens Observa- 
tory the magnitude is Ms(ATH)=ML+0. 5=6.0; we clar- 
ify that predictions are issued only when the expected 
magnitude is comparable to 6.0 units or larger). Note 
that this earthquake occurred after the initial submis- 
sion of the present paper (and its time of occurrence was 
approached 54 by means of the procedure developed in 
Ref^). 

Here, we analyze, as an example, the long duration 
SES activity that lasted from February 29 until March 
2, 2008(FigJ3][c)). The time-series of this electrical dis- 
turbance, which is not of obvious dichotomous nature, 
is depicted in the upper channel (labelled "a") of FigJU 
The signal, comprising a number of pulses, is superim- 
posed on a background which exhibits frequent small MT 
variations. The latter are simultaneously recorded at all 
measuring sites, in contrast to the SES activities that 
solely appear at a restricted number of sites depending 
on the epicentral region of the future earthquake 3 ^.. This 
difference provides a key for their distinction. In order to 
separate the background, we proceed into the following 
steps: First, we make use of another measuring dipole of 
the same station, i.e., the channel labelled "b" in Figj4j 
which has not recorded the signal but does show the MT 
pseudo-sinusoidal variations. Second, since the sampling 
rate of the measurements is 1 sample/ s, we now compute 
the increments every Is of the two time-series of channels 
"a" and "b" . Assuming the "Is increments" of channel 
"a" lying along the x-axis and those of "b" along the y- 
axis, we plot in the middle panel of Figjl] (labeled "c") 
the angle of the resulting vector with dots. When a non- 
MT variation (e.g., a dichotomous pulse) appears (disap- 
pears) in channel "a" , this angle in "c" abruptly changes 
to 0° (±180°). Thus, the dots in panel "c" mark such 
changes. In other words, an increased density of dots 
(dark regions) around 0° or ±180° marks the occurrence 
of these pulses, on which we should focus. To this end, 
we plot in channel "d" of FigH] the residual of a linear 
least squares fit of channel "a" with respect to channel 
"b" . Comparing channel "d" with channel "a" , we notice 
a significant reduction of the MT background but not of 
the signal. The small variations of the MT background 
that still remain in "d" , are now marked by the light blue 
line, which when removed result in channel "e" of FigJU 
Hence, channel "e" solely contains the electric field vari- 
ations that precede rupture. This channel provides the 
time-series which should now be analyzed. 

The DFA analysis (in conventional time) of the time- 
series of channel "e" of FigJU is shown in FigJSJ It 
reveals an almost linear log F(n) vs logn plot (where 
n = f exp At) with an exponent a « 1 practically over all 
scales available (approximately five orders of magnitude) . 
Note, that this value of the exponent remains the same 
irrespective if we apply DFA-1, DFA-2 or DFA-3. This 
result is in agreement with the results obtaine d 25 ! 26 ' 40 ! 41 
for SES activities of shorter duration. 

For the classification of this signal, i.e., to distinguish 



whether it is a true SES activity or a man-made electric 
signal, we now proceed to its analysis in natural time. In 
order to define (xk, Qk) the individual pulses of the sig- 
nal depicted in channel "e" of Fig(?]have to be identified. 
A pulse starts, of course, when the amplitude exceeds a 
given threshold and ends when the amplitude falls below 
it. Moreover, since the signal is not obviously dichoto- 
mous, instead of finding the duration of each pulse, one 
should sum the "instantaneous power" during the pulse 
duration in order to find Qk, as already mentioned. To 
this end, we plot in Fig(B] the histogram of the instan- 
taneous power P of channel "e" of Fig01 computed by 
squaring its amplitude. An inspection of this figure re- 
veals a bimodal feature which signifies the periods of in- 
activity (P< 500^V 2 Hz) and activity (P> 500AiV 2 Hz) 
in channel "e" of FigJ3J In order to find Qk, we focus 
on the periods of activity and select the power threshold 
Pthres around the second peak of the histogram in Fig|5] 
Let us consider, for example, the case of Pt/ lres =1400 
^V 2 Hz. In FigJTJa), we depict the instantaneous power 
P of the signal in channel "e" of FigJU Starting from 
the beginning of the signal, we compare P with Pthres 
and when P exceeds Pthres we start summing up the P 
values until P falls below Pthres for the first time, k = 1. 
The resulting sum corresponds to Q\. This procedure is 
repeated until P falls below Pt/ lres for the second time, 
k = 2, and the new sum represents Qi, etc. This leads 
to the natural time representation depicted in Fig [7Kb). 
The result depends, of course, on the proper selection of 
Pthres ■ The latter can be verified by checking whether a 
small change of Pthres around the second peak of the his- 
togram, leads to a natural time representation resulting 
in approximately the same values of the parameters n\, S 
and S- . By randomly selecting Pthres in the range 500 to 
2000 /iV 2 Hz, we obtain that the number of pulses in chan- 
nel "e" of Figgjis N = 1100±500 with /t x = 0.070±0.007, 
S = 0.082 ± 0.012 and S- = 0.078 ± 0.006. When P thres 
ranges between 1000 to 1500 /xV 2 Hz, the corresponding 
values are N = 1200 ± 200 with ki = 0.068 ± 0.003, 
S = 0.080 ±0.005 and S- = 0.074 ±0.003. Thus, we ob- 
serve that irrespective of the Pt/ lres value chosen, the pa- 
rameters Ki, S and S- obey the conditions (J5J) and © for 
the classification of this signal as SES activity. The same 
holds for a non-dichotomous signal recorded on March 28, 
2009 at Keratea station located close to Athens (FigJSJ). 
To approach the occurrence time of the impending event, 
the procedure developed in Refi22 has been employed for 
the seismicity within the area Ngy yE^ g. The results, 
when considering the seismicity until early in the morn- 
ing on June 19, 2009, are shown in FigO They show 
that a local maximum of the probability 32 Prober ) of 
the seismicity is observed at k\ = 0.070 upon the occur- 
rence of the following events: (a) an Ml=3.2 event at 
12:42 UT on June 18 with epicenter at 38. 7N 23.0E and 
an M L =3.2 event at 03:44 UT on June 19 at 37.9N 23.0E 
if we take a magnitude threshold M t /j res =3.0 (b) when 
taking Mthres—S.l, the maximum occurs again upon the 
occurrence of the aforementioned event at 12:42 on June 
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18 (c) when increasing the threshold at M t / lres =3.2, the 
maximum occurs somewhat earlier, i.e., upon the oc- 
currence of the Mi=3.5 event at 08:26 UT on June 18 
at 38. 7N 23. 0E. Let us now consider a larger area, i.e., 
NlflesEffioj an d repeat the same calculation. The re- 
sults are shown in FigfTUl They reveal again the occur- 
rence of a local maximum of Prob(«i) at K\ — 0.070 upon 
the occurrence of the above mentioned Ml=3.2 event at 
12:42 UT on June 18 when considering either M t ^ res =3.2 
(FigJIOTc)) or M^ res =3.1 (FigflOTb)) (for the latter case 
a maximum also occurs at the last event shown in the 
figure with M L =3.2 at 03:44 UT on June 19). The 
maximum appears upon the occurrence of a Mi=3.0 at 
23:48 UT on June 18 when decreasing the threshold at 
M thres =3.0. 

In summary, the results exhibit both magnitude- and 
spatial- invariance, which is characteristic when ap- 
proaching the critical point. 

Note added on July 21, 2009. Repeating the aforemen- 
tioned calculation for the area Nl^ggEjf^o, we find the 
results depicted in FigfTlTa) for M f /j res =2.8. An inspec- 
tion of this figure reveals a local maximum of Prob(/ii) 
at Ki = 0.070 upon the occurrence of a Mi=3.0 at 21:10 
UT on July 19, 2009 with an epicenter at 38.1N 22. 7E 
(the same result is found for Mt; 1I , es =2.9.) Interestingly, 
this local maximum turns to be the prominent one when 
studying the seismicity in the same area (for the same 
threshold M tferes =2.8) since June 21, 2009, see FigfJUb); 
the same holds for M t h, res =3.0. Furthermore, the lat- 
ter study was made for the smaller area NgyfE^I for 
Mt/j res = 2.6, 2.7 and 3.0 and showed a local maximum 
of Prob(Ki) at Ki — 0.070 upon the occurrence of the 
same event (i.e., at 21:10 UT on July 19, 2009). The in- 
vestigation still continues in order to clarify whether the 
latter local maximum will eventually turn to become a 
prominent one. 

Note added on September 5, 2009. At 8:17 UT on 
September 2, 2009, a Mi=4.3 earthquake occurred al- 
most at the center of the predicted area N^yyEi^'g. The 
continuation of the investigation of the previous para- 
graph in the smaller area NgyyE^'l showed a prominent 
maximum of Prob(Ki) at n\ = 0.070 (see FigfT2] (a)) 
upon the occurrence of a Mi=3.3 event at 09:35 UT on 
September 2, 2009, with an epicenter at 38.1N 23.3E. 
This was found to hold for the aforementioned larger area 
(see FigfT2](b)) as well, thus revealing that the criticality 
condition -as developed in Ref. l32l - is obeyed. 

Note added on September 19, 2009. The previous Note 
was followed by three earthquakes of magnitude 4.4, 4.7 
and 4.8 that occurred at 21:15 UT September 9, 06:04 
UT September 10 and 07:12 UT September 16, 2009 with 
epicenters at 38. 7N 22.8E, 38.3N 24.1E and 39.0N 22. 2E, 
respectively. The continuation of the natural time inves- 
tigation of the seismicity for M t /j res =3.0 in the smaller 
area NgyyE^f 'g revealed a prominent maximum Prob(Ki) 
at ki — 0.070 (see FigfT3] (a)) upon the occurrence of a 
M L =3.0 event at 01:55 UT on September 18, 2009, with 
epicenter at 38.7N22.8E. Since the same was found to be 



true for the aforementioned larger area as well (see FigfT3l 
(b)), the investigation is now extended to additional mag- 
nitude thresholds in order to clarify whether (beyond 
spatial and) magnitude threshold invariance holds, thus 
identifying that the system approaches the critical point. 

Note added on October 13, 2009. The previous note 
was actually followed by an earthquake of magnitude 
4.6 that occurred at 22:27 UT on October 2, 2009 with 
an epicenter at 38. 3N 22. 9E. The natural time investi- 
gation of the seismicity still continues with the follow- 
ing results. First, in the smaller area Ngy-yEjf-g: For 
M t hres=2.7, the criticality condition started to be ap- 
proached since October 7 and the most prominent max- 
imum of Prob(fti) at n\ = 0.070 was observed upon the 
occurrence of the Ml=2.7 event at 01:56 UT on October 
11, 2009. The same behavior was found for M t /j res =2.8 
(but the last event for this case occurred at 08:15UT on 
October 10, 2009 at 38. 2N 23.1E). Second, in the larger 
area N^^E^;^: For both thresholds M thres =2.7 and 
M t hres=2-8, the criticality condition was also started to 
be approached since 20:18 UT on October 7, 2009. The 
most prominent maximum of Prob(fti) at k\ = 0.070 for 
M t f l res='2-7 have been observed upon the occurrence of 
the aforementioned M^=2.7 event at 01:56 UT on Oc- 
tober 11, 2009, while for M t / lres =2.8 it was observed at 
4:51 UT on October 9, 2009, upon the occurrence of the 
M L =2.9 event at 38. IE 23.2E. Interestingly, the results 
for the larger area are almost identical if the investigation 
starts at 15:00 UT on July 31, 2009. 

Note added on November 14, 2009. Meanwhile, two 
SES activities were recorded at PAT on October 24 and 
November 11, 2009, see FigfMl Hence, the ongoing study 
(Fig iT5|) of the previous Note should be accompanied by 
a complementary study, focused on the determination of 
the occurrence time of the impending mainshock, that 
also considers the seismicity in the western part (i.e., 
NjjyjEfgJ) of the PAT selectivity map (see area A in 
Fig. 8 of Ref.—). Such a complementary study, when 
taking into account the seismic data until early in the 
morning on November 13, 2009, leads to the results de- 
picted in Fig ll6l for Af t ^ res =3.1. A local maximum of 
Prob(Ki) at Ki = 0.070 was observed (see the rightmost 
arrow in Fig lTB"]) upon the occurrence of Ml=3.2 event 
at 18:55 UT on November 2, 2009 (that was followed by 
a M L =5.6 EQ at 05:25 UT on November 3, 2009, with 
an epicenter at 37.4°N 20.4°E). In addition, the most 
prominent maximum in Fig |16l of Prob(Ki) at k\ — 0.070 
is observed later, i.e., upon the occurrence of a Mi=3.1 
event at 01:42 UT on November 13, 2009 (see the left- 
most arrow in Fig |16l) . This, which has been also checked 
for M t hres—^-0, provides evidence that the critical point 
is approached. Calculations are now carried out to fur- 
ther investigate the magnitude threshold invariance of 
this behavior. 

Note added on November 27, 2009. Actually, after the 
previous Note the following two earthquakes occurred at 
18:30 UT on November 16, 2009 and at 04:29 UT on 
November 19, 2009 with magnitudes (M S (ATH)) 4.3 and 
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4.4 and epicenters at 38.4°N 22.0°E and 38.2°N 22.7°E, 
respectively, that both lie within the PAT selectivity 
map, i.e, NpjEfgJ. The continuation of the two studies 
mentioned in the previous Note, i.e., in the PAT selectiv- 
ity map as well as in the area N^-ggE^f -J,, both show that 
the most prominent maximum of Prob(Ki) at Ki = 0.070 
was observed upon the occurrence of the M^=3.2 event 
at 19:31 UT on November 25, 2009, thus indicating that 
the system is near the critical point. 

Note added on February 9, 2010. Two earthquakes 
with magnitudes M S (ATH)(= Ml + 0.5)=5.7 and 5.6 oc- 
curred at 38.4°N 22.0°E (and hence inside the expected 
area) at 15:56 UT on January 18 and at 00:46 UT on Jan- 
uary 22, 2010 (after observing a maximum of Prob(Ki) 
versus n\ at ki=0.070 for M t ^ res =3.1 upon the occur- 
rence of the Mi=3.2 event at 15:56 UT on January 15, 
2010). These two magnitudes are consistent with the con- 
dition M S (ATH)>6.0 (under which SES information are 
publicized) if a reasonable^ experimental error 3<r=0.7 
is accepted. 

However, since the amplitudes of all the SES activ- 
ities recorded could warrant an impending mainshock 
with magnitude exceeding M S (ATH)=6.0, and in order 
to determine the occurrence time of such a mainshock, if 
any, the natural time analysis of seismicity still contin- 
ues. This analysis starting at 04:00 UT on December 27, 
2009, revealed that in the area N||;gE|f;| (almost 200km 
x 100km) the criticality condition (i.e., Prob(Ki) ver- 
sus K\ maximizes at ki=0.070) for M t / lre s=3.5 is obeyed 
during the last three events of FigflTl the latest of which 
(that exhibits the strongest maximum) occurred at 10:18 
UT on February 7, 2010 at 38.6°N 23.7°E with Mi, =4.3. 
This has been checked to be valid for broader areas as 
well (spatial invariance), but has not yet been found to 
be confirmed (until early in the morning on February 
8, 2010) for lower magnitude thresholds. The investiga- 
tion of the latter is still in progress, because the following 
might happen: Due to coarse-graining the criticality con- 
dition might have been obeyed at a larger time-window 
for large magnitude thresholds, while for smaller thresh- 
olds the time- window may become smaller thus achieving 
a better accuracy in the determination of the occurrence 
time of the impending earthquake. 

IV. DISCUSSION 

In general, electric field variations are interconnected 
with the magnetic field ones through Maxwell equations. 
Thus, it is intuitively expected that when the former ex- 
hibit long range correlations the same should hold for the 
latter. This expectation is consistent with the present 
findings which show that, at long time-lags, the original 
time-series of both electric- and magnetic field variations 



preceding rupture exhibit DFA-exponents close to unity. 

The above can be verified when data of both electric- 
and magnetic-field variations are simultaneously avail- 
able. This was the case of the data presented in FigfT] 
In several occasions, however, as mentioned in Section I, 
only magnetic field data exist, in view of the fact that it is 
easier to conduct magnetic field measurements compared 
to the electric field ones. When using coil magnetome- 
ters, the magnetic field variations have the form of series 
of 'spikes'. Whenever the amplitude of these 'spikes' sig- 
nificantly exceed the pseudo-sinusoidal variations of the 
MT background, as in the case of Bew in Fig[TJ a di- 
rect application of DFA (see Figj2]) elucidates the long 
range correlations in the magnetic field variations pre- 
ceding rupture. On the other hand, when considerable 
pseudo-sinusoidal MT variations are superimposed, a di- 
rect application of DFA is not advisable. One must first 
subtract the MT variations (following a procedure sim- 
ilar to that used in the electric field data in Fig0]) and 
then apply DFA. 

The preceding paragraph refers to the analysis of the 
signal in conventional time. As already shown in Ref. 41 , 
natural time analysis, has the privilege that allows the 
distinction between true SES activities and man made 
signals. This type of analysis, however, demands the 
knowledge of the energy released during each consecu- 
tive event. The determination of this energy is easier 
to conduct in the case of electric field variations (this 
is so, because coil magnetometers, as mentioned in Sec- 
tion II, act as dB/dt detectors). When these variations 
are of clear dichotomous nature, the energy release is 
proportional to the duration of each puls o 40 ' 41 . On the 
other hand, in absence of an obvious dichotomous nature 
an analysis of the 'instantaneous power' similar to that 
presented in the last paragraph of Section III should be 
carried out to determine the parameters K\, S and S_ in 
natural time. 



V. CONCLUSIONS 

First, DFA was used as a scaling analysis method to 
investigate long-range correlations in the original time 
series of the magnetic field variations that precede rup- 
ture and have the form of 'spikes' of alternating sign. We 
find a scaling exponent a close to 0.9 for time-lags larger 
than the average time interval (T) between consecutive 
'spikes', while at shorter time-lags the exponent is close 
to 0.5 thus corresponding to uncorrelated behavior. 

Second, using electric field data of long duration SES 
activities (i.e., from several hours to a couple of days) 
recorded during the last few years, DFA reveals a scale 
invariant feature with an exponent awl, over all scales 
available (approximately five orders of magnitude). 
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FIG. 1: (a): Variations of the electric- (upper three channels) 
and magnetic- (lower two channels) field recorded on April 
18, 1995 (see the text). (b):A lOmin excerpt of (a). The two 
dashed lines in (a) show the excerpt depicted in (b). 
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FIG. 2: (color online)The DFA for the B E w channel of 
Fig. 1(a). Logarithm base 10 is used. 
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FIG. 3: (color online) Long duration SES activities during 
the last few years (see the text). 
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FIG. 4: (color online) The long duration SES activity from 
February 29 to March 2, 2008. Channel "a":original time 
series, channel "b": recordings at a measuring dipole which 
did not record the SES activity, but does show MT variations, 
"c" :the angle of the resulting vector upon assuming that the 
"Is increments" of channel "a" lie along the x-axis and those 
of channel "b" along the y-axis. Channel "d" : the residual of a 
linear least squares fit of channel "a" with respect to channel 
"b", channel "e":the same as "d" but after eliminating the 
slight variations of the background. For the sake of clarity, 
channels "a" , "b" and "d" have been shifted vertically. 
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FIG. 5: (color online) The DFA-Z (1=1,2 and 3) for the lower 
channel (i.e., the one labelled "e") of Fig|4] Logarithm base 
10 is used. 
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FIG. 6: Histogram of the "instantaneous power" P, i.e., the 
squared amplitude of the signal depicted in channel "e" of 
FigH 
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FIG. 7: (color online) (a) The "instantaneous power" P of the 
signal depicted in channel "e" of Fig|4] (computed by squaring 
its amplitude). The solid line parallel to the x-axis marks an 
example of a threshold Pih res (=1400 /iV 2 Hz) chosen, (b) 
The resulting representation of the signal depicted in channel 
"e" of Fig|4]in natural time, when considering Pth res = 1400 
MV 2 Hz. 




200 400 600 800 1000 1200 1400 

Time(s) 



FIG. 8: The most recent electric signal analyzed (see the 
text). 
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FIG. 9: The probability Prob(Ki) of the seismicity within 
the area Nf^f E||;g subsequent to the SES activity recorded 
at KER on March 28, 2009, when considering the following 
magnitude thresholds: M t(lres =3.0 (a), 3.1 (b) and 3.2 (c). 
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FIG. 10: The same as FigO but for the area N|?;^Elfig in 
order to check the spatial invariance of the results (see the 
text). 
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FIG. 11: The same as FigU but for the area Np;^El|;^ by 
considering the seismicity since (a) March 28, 2009 and (b) 
June 21, 2009 (see the text). 
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FIG. 12: The probability Prob(/-ci) vs /ti upon considering 
the seismicity (M(/ tr . es =3.1) until 9:35 UT on September 2, 
2009: Panel (a) corresponds to the area Ni?:?E||i, whereas 
panel (b) to the larger area N^gsEl^o- 
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FIG. 13: The probability Prob(Ki) vs m upon considering 
the seismicity (Mthres~3.0) until 1:55 UT on September 18, 
2009: Panel (a) corresponds to the area Ni?:?Eiy, whereas 
panel (b) to the larger area N|f:ggE|i^. 
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FIG. 14: The SES activities at PAT on 24 October 2009 (up- 
per) and 11 November 2009 (lower), see the text. 
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FIG. 15: Earthquakes with magnitude Ms(ATH)>4.5 and 
epicenters within (or very close to) the area(s) specified in 
advance that followed the consecutive Notes added to this 



paper. 




FIG. 16: The probability Prob(fti) of the seismicity within 
the area Nf^jjEjJJ subsequent to the SES activity recorded 
at PAT on October 24, 2009, when considering M t h res =3.1. 
The seismic data until early in the morning on November 13, 
2009, have been used. 




FIG. 17: The probability Prob(/-ci) of the seismicity, for 
M t hres=3.5, within the area NH^E^iis when considering the 
seismic data until early in the morning on February 8, 2010. 



